10 PRINT "EXAMPLE 1.4-2 ENERGY LOS5 IN SLEEVE": PRINT
20 CLEAR : DEFINT I, N: DIM F(36), X(36)
30 DEF FNX (TH) = R * SQR(1! - (RR * SIN(TH) / R) ^ 2) - RR * COS(TH)
40 DEF FNV (TH) = OM * RR * SIN(TH) * (1! - RR * COS(TH) / SQR(R ^ 2 - (RR * SIN(TH)) ^ 2))
50 READ R, RR, D, L, MU, DELTA, RPM, N, PI
60 DATA 2., .5, .8, 8., 1E-4, .001, 1200., 36, 3.1416
70 PRINT "R,RR,D,L="; R; RR; D; L: PRINT "MU,DELTA,RPM,N,PI="; MU; DELTA; RPM; N; PI
80 OM = 2! * PI * RPM / 60!: PERIOD = 2! * PI / OM: DT = PERIOD / (2! * N)
90 C1 = MU * PI * D * L / (12! * DELTA): W = 0!: 'FORCE=F=C1*V
100 FOR I = 0 TO N: TH1 = I * OM * DT: X(I) = FNX(TH1): F(I) = C1 * FNV(TH1): NEXT I
110 FOR I = 1 TO N: 'TRAPEZODAL RULE INTEGRAT10N FOR HALF PERIOD
120 W = W + .5 * (F(I) + F(I - 1)) * (X(I) - X(I - 1)): NEXT I
130 PRINT : POWER = W / (.5 * PERIOD)
140 PRINT "POWER="; POWER; " FT-LB/S"